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Introduction. The solution to the problem on the centroidal 
motion control synthesis (guidance problem) of an unmanned 
aerial vehicle (UAV) with long-range capabilities in the boost 
phase is considered. Control condition requires optimum fuel 
consumption. The principle of dynamic programming consid- 
ering the restrictions to the vector modulus of the thrust output 
is used to solve the problem. The implementation of terminal 
guidance requires the formation of control as a function of the 
object state at the end of the ascent phase. The attainment of 
these boundary conditions determines the further transition to 
the ballistic flight phase. 

Materials and Methods. Bellman’s principle of dynamic pro- 
gramming is the most reasonable from the point of view of the 
implementability of the computationally efficient on-board 
algorithms and the solution to the problems in the form of 
synthesis. With natural scarcity of thrust and energy resources 
on board, this principle enables to obtain solutions free from 
the switching functions. In this case, the optimal control is a 
smooth function (without derivative discontinuity) of the cur- 
rent and final parameters of the UAV. 

Research Results. A new algorithmic method for the synthesis 
of terminal motion control is developed. Its difference is that 
the UAV movement control in the ascent phase is formed by 
the function of the motion actual and terminal parameters. 
This ensures movement along an energetically optimal 
trajectory into the given region of space. The problem solution 


results enable to build closed terminal guidance algorithms for 


” The research is done within the frame of the independent R&D. 


Beedenue. CratTba MOcBAIIeHa pellleHHIO 3aya4H CHHTe3a 
ylipaBJIeHHa JIBYXKeHHeM IleHTpa Macc (3aqa4a HaBeyleHHs) 
OecHHIOTHOrO JIeTaTesbHOrO amnapata (bJIA) c OonpmoL1 
JaJIbHOCTbIO MOJIeTa Ha pa3srOHHOM yyacTKe. YcuIOBHe yupaB- 
JICHHA: ONTHMAJIbHbIM pacxoy TomIMBa. Jia pelieHua 3a7aqn 
HCHOUIb3yYeTCA TIPHHUMI THHAaMMYecKOrO MporpamMMupoBaHuA 
C y4eTOM OrpaHHyeHHt Ha MOJlyJIb BEKTOpa TAM JBUTaTessA. 
Peamn3allua TEPMHHAJIbHOrO HaBeyeHuA TpeOyeT CPOpMupoBa- 
HUA yiIpaBleHud Kak (yHKUMM COCTOAHUA OObEKTAa B KOHIIC 
yuacTka BbIBeyeHua. JlocTuxKeHHe 3THX TpaHH4HbIX ycuOBUi 
olpeyeiaeT TaibHeiiiMu mepexoy kK OammMcTHYecKON da3ze 
moyleTa. 

Mamepuaazoi u memoooi. Wpunuun uHamuyeckoro mporpam- 
MupoBaHna beswimMaHa ABJIAeTCA HaMOOIIee pallMOHaJIbHbIM C 
TOUKH 3PeCHHA peallMu3yeMOcTH IPeKTHBHBIX B BbIYHCIIMTEIb- 
HOM OTHOMICHHH OOPTOBHIX aJITOPHTMOB UH pellleHHA 3aadu B 
(bopMe cuHtTe3a. IIpu ecTrecTBeHHOM OrpaHHyeHHOCTH BeJIM- 
4HHbI TATH HW 9HEPreTH4YeCKMX PecypcoB Ha OoOpTy MaHHbIi 
TIpHHUMM MO3BOJAeT MOJYYINTb pellleHua, He coOseprxkaHe 
(@YHKUMU MepekHOueHHA. ONTHMasIbHOe yiIpaBlIeHve B 3TOM 
ciyyae ABJIAeTCA TaKOU PyHKUMet (Oe3 pa3spbiBa IpOu3BO- 
HOM) TeKYUIMX UW KOHCUHBIX MapamMeTpos BJIA. 

uccIeOOBAHUA. 

clocoo 


Pezylbmambl Pa3pa0oTaH HOBBIM 


asIrOopHTMuYecKHi CHHTe3a TePpMMHasIbHOrO 
yipaBsieHuad BMyKeHuA. Ero oTIM4Me B TOM, 4TO yiIpaBsIeHHe 
mpwxkeHnueM BbJIA Ha pa3roHHOM y4acTKe TpaeKTopHu 
(bopmupyetca PyHKIMeH TeKyUIMX MH KOHCYHbIX MapaMeTpoB 
mBwxKeHua. Takum o0pa30M oOecneuMBaeTCA JIBYHXKeHHe TO 
9HepreTHyeCKH OMTHMaJIbHOM TpaeKTOpuH B  3aaHHyto 
oOmacTb wWpoctpaHctBa. Pe3yiIbTaTbl pelleHua 3aa4u 
TIO3BOJIAIOT CTPOUTb 3AMKHYTbIC aJITOPUTMbI TePMHHaJIbHOrO 
HaBeyIeHHA JIA pa3rOHHOrO yyacTKa TpaekTopuu BJIA c 


OObIION WaIbHOCTbIO MoseTa. Takve asIrOpHTMbI OOaTaroT 
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the boost phase of the UAV trajectory with long-range 
capabilities. Such algorithms have good convergence and 
injection accuracy due to the prediction of parameters during 
the flight at a shorter time interval. 

Discussion and Conclusions. The most preferred is_ the 
principle of dynamic programming. It should be used when 


solving the problem on the centroidal motion control synthesis 


XOPOMIeH CXOJMMOCTbIO MH TOUHOCTHIO BbIBeCJICHHA 3a CUCT 
IIPOrHO3HpOBaHHA TapaMeTpoB B porecce oseTa Ha 
COKpalljarollleMca HHTepBasie BPeMeHH. 

O6écyacdenue u 3akmoueHua. Handomee lpexqnNOuTUTeIbHbIM 
IIpecTaBsAeTca IIpHHIMH MHAMUYeCKOrO 
porpamMmMupospanna. MMeHHo ero cileyeT MCHOJIb30BaTb TIPU 
pelleHHu 3aja4dn CHHTe3a ONTMMAIIbHOrO MO pacxoyy TOMIMBa 


ylipaBJIeHHa JIBHXKeHHeM IleHTpa Macc (3aqa4a HaBejleHHs) 


(guidance problem) of the UAV with long-range capabilities BJIA c OoMbUIO WaIbHOCTbIO HoseTa Ha pa3srOHHOM y4acTKe. 


in the boost phase. 


Keywords: unmanned aerial vehicle (UAV), terminal guid- KuroueBbie CcJI0Ba: OeCIIMJIOTHBIM JeTaTeIbHbIM arilapat 
ance, direction cosines, pitching angle, angle of attack, bound- 


ary conditions, boost phase, ballistic flight phase. 
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Introduction. Currently, the capabilities and scope of applicability of the unmanned aerial vehicles (UAV) 
have increased significantly. This is primarily due to the UAV flying range extension. Vehicles with a ballistic phase of 
flight should be provided with the control in various areas, including the UAV positioning. For the control design, it is 
reasonable to apply the principle of terminal guidance. The control should be developed as a function of the terminal 
motion variables, and not as a temporal function. 

A great many publications are devoted to this problem solution, but the task described above holds relevance. 
In particular, the development of highly computationally efficient algorithmic methods of terminal guidance 1s of inter- 
est. At the same time, the features universal for various types of launch vehicles should be considered. They have adap- 
tive features and in a certain sense meet the optimality requirements. 

Considering the guidance task, such a control, which uses the minimum amount of fuel, is valuable. The syn- 
thesis of optimal control is based on the use of the Bellman dynamic programming method [1]. 

Materials and Methods. The solution to the problem of the synthesis of optimal control of the UAV motion in 
the boost phase of power-on flight has been studied in many papers. However, the optimal control solution obtained in 
most cases is reduced to the implementation of a time or parametric program. We will solve the problem of finding the 
optimal control in the UAV ascent phase in the following formulation. The UAV motion parameters are known as Xo, 
Yo, Zy) coordinates of the current path point, obtained through solving the problem of navigation. X,, Y,, Z, parameters of 
the final point satisfy the boundary state that fixes the transition to the ballistic phase of flight determined by the hyper- 
surface in the phase space. The condition of S; (X,., Y,, Z,.) = 0 1s satisfied by a whole set of finite parameters. It is re- 
quired to synthesize optimal control in the task of UAV guidance, which ensures its transfer from the initial state to the 
hypersurface of the final conditions. 

The optimality criterion is the amount of fuel consumed in the boost phase: 


m(t) = ri(t)dt (1) 


where m(t) 1s fuel mass flow rate. 


The following system of differential equations is adopted as a mathematical model of the centroidal motion: 
R®=Ve), (2) 
V(tV=WOEO+ 87), 
where R(f) is radius-vector, Z(7) is terrestrial attraction vector, W(t) is module of control acceleration vector. 
The module of the control acceleration vector (W(t) ) is a specified temporal function, and it is determined by the 


UAV engine performance. Unknown is the unit vector of control acceleration — E(t). When solving the optimization 


problem, it will determine the required properties of the UAV guided motion in the boost phase. There is no need in the 
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thrust vector control (in a sense, this is apparent acceleration amount — W(t) ). A more rational approach is to maximize 


this value for the engine booster. This ensures the application of the solutions obtained for the case of using solid fuel 
engines. The UAV engine performance is stable enough, and at m(t) = const constant fuel-flow rate, the optimal con- 


trol is determined by the combustion duration. In this case, the functional (1) will be a function of the upper limit of 
integration. Hence, the problem of minimizing the fuel amount turns into an equivalent task of minimizing the flight 
time, and speed-of-response will be the optimality criterion. The task of the synthesis is to find the orientation of the 


thrust vector of the UAV engine, which is determined by the direction cosines of the thrust vector ( P(t) ) as a function 


of the current parameters and the final state. 

We admit two assumptions. The first is as follows. Since a large part of the boost phase lies outside the dense 
atmosphere, we will not consider the angle rate of evolution of the thrust vector in space (longitudinal axis of the UAV) 
as constraint. The second assumption is the following. Assume that the UAV motion is passing in a predetermined 
plane. 

The UAV final state for the transition to the ballistic phase of flight is fixed by satisfying the following bound- 
ary condition [2]: 

Si = (Vy, =x )[V (y, -y,)-V,, ie —%, )|- 
XX TV, (3) 


U1 y k k 


Here, x,,),,X,,¥, are, respectively, the coordinates of the initial point and the point of the starting of the UAV opera- 


tion on the final path segment. The current value of the boundary condition (S,,(¢,) ) is a measure of non-compliance 


with the condition (3). A mathematical notation for this condition corresponds to the hypersurface, which is a smooth 
function of phase coordinates and describes the entire family of possible UAV ascent trajectories [2]. 

The boundary conditions are specified for the central field of attraction. For this case, the equations of the UAV 
motion will have the following form: 


x, =%as 
Tt 
{= XP COSC, 
"x =x, 4) 
: , 
Tt 
Xx, =-—x,+—Pcosa,. 
3 3 2 
r m 


Here, =, = f-M, 1s constant of the central field of the earth’s attraction equal to the product of the gravitational con- 


stant (f) and the mass of the Earth (M,); r=(x; +x, )'", x, =x, x,=V,,x,=y,x,=V, ; P is thrust vector value; X4 is aero- 


x9 "3 
dynamic drag force; Y, is aerodynamic lift. 
Since the UAV is fitted up with a steerable thruster, its direction cosines will determine its pointing. In this case, 
the handling constraint will be determined by the ratio: 


|PO|={[PO cosa, (OF +[PO cose, One <|P°(o). (5) 


The handling constraint (5) is “hypersphere restriction”. This implies the solution in which optimal control is 
not piecewise constant, without switching. In this case, the time optimal ascent trajectory in the phase space has no “an- 
gles” and no discontinuities of the derivative. The time optimal control will be a nonlinear, continuous function of the 
boundary conditions (3) [3]. 

The synthesis task is formulated as follows. The control object from an arbitrary current state, taken as the ini- 
tial one and determined by the current value (Sj;,) at the time (f), 1s transferred to the hypersurface of the boundary con- 
dition S|, =0 at the time (¢,). This takes into account the control acceleration amount (engine thrust module). It is re- 
quired to find the optimal control in the form of synthesis, which provides such a transfer in the shortest possible time. 

The plant state at the final instance satisfies the boundary condition S,[X(t,)]=0 and determines the moment 


of transition to the ballistic phase of flight. 
In accordance with R. Bellman’s dynamic programming principle [4], the necessary and sufficient condition 
for optimality of the formulated problem will be the ratio: 


nip] 1.0) =-], (6) 


ueU i=l 


Considering the mathematical model of the object specified by the system (4), the dynamic programming equa- 
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tion [4] with the optimality criterion, is determined by the relation: 


0 0 0 0 
np a [Fe - tp |S +x ca -(3s-22)2 |-+. (7) 


a 4 
Oe Ox, r m Ox, Ox, r m Ox, 














where P, = Pcosa,,P, = Pcosa,. 


Due to the methodology of Bellman’s dynamic programming, the minimization can be carried out through the 
application of the Schwartz inequality [3] to the relation (5). This will significantly simplify the solution to the optim1- 
zation problem. Then, the expression (7) with regard to inequality (5) will be determined by the relation: 














2 2 1/2 
0 0 0 0 
gi oe | | (8) 
wU | m ~ Ox, Mm” OX, m OX, m OX, 
In this case, the Hamilton-Jacobi equation will be as follows: 
OT’ 1%, OT” 








Tt 
0 0 
x ee, ye (LX 
Ox, I Ox, 1 


: ae ; 9 
1 oT 1 oT ) 
x} | — +| — =I, 
m Ox, m OX, 


The equation (9) is solved considering the given boundary condition: 























(X15 Xp> Xo %,) ES - (10) 
The calculation of the partial derivatives of S for phase variables (x,,x,,x,,x,) gives the following dependenc- 
es: 
os” 
— (-2x,x,x, nee Wee 90 ee ee 
1 
os* ; 
Fs = (2x; Ky 2X Ae 2K Xi ON Xa, ) 
<a (11) 
os ; ; 
Ay = wo KH 2L Me Ht HW, ), 
3 
os* ; 
rs = (2x; Ki ALN LK KG PAH, PX NX ). 
4 
Through substituting the relations (11) into (9) and performing simple transformations, we obtain the depend- 
ence: 
P 0 
{= @, SR AZ N i, oe 
r m : 
x{ (2x x, LU NN LK Ns NN +X,%5%,, ) + (12) 
or: 
+(2x7%, SZ Oh a, ) lian — —] 
1 
with the boundary conditions 
T°’ (S,) =0, at (x,, x,,%,,x,) €5,,. (13) 
The expression in curly brackets in (12), denoted by S’(x,), can be written in a compact form: 
oT” 
S’(x.)- =-—]. 14 
(x; ) as (14) 


The equation (14) with boundary conditions (13) can be solved in various ways, for example, by the method of 
characteristics [5]. However, it is more rational to use the following method. 
The expression (11) defines the optimal control structure [6]: 


ae (x) = —PRe”" (x), 


oe - (15) 
We (x) — WE" (x), 
where E(x) is unit thrust vector. 
The ratio (15) in the scalar form: 
P” = Pcosa,, 
P™ = Pcosa,, ~ 
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ar® 
Ox, 
COSQ, = rT, Ung 2? 
OT” oT” 
|| Sy 
Ox, Ox, 
ax, 
COSQ, = 


2 2 1/2 * 
OT” OT” 
—— + —__ 
Ox, Ox, 


The equation (9) and relations (17) can be put in a more convenient form for the following reasons. According 
to the formation of a set of perturbed trajectories, through varying the control, it is possible to construct a hypersurface 
with equal ascent time, that is, an isochronous surface, in the phase space. Indeed, for each phase trajectory point, we 
will calculate the optimality criterion value: J(t,), J(t),...,/(t). 

Hence, we obtain a set of trajectories for each fe [t,t]. Owing to the continuity of the dependence of x(¢) and J 
on the variable control, a set of trajectories forms a surface in the phase space (X). This boundary surface formed by the 
set of vectors x[t;,/(t;)], 1s convex and smooth. For transition conditions to the ballistic flight phase, the isochronous 
surface has a tangency point with the hypersurface of the boundary conditions. 

Under the qualitative implementation of the optimal control, the distance in the phase space between the hyper- 
surface of the boundary conditions S;=0 and the isochronous surface T (x, R,, t) = 0 will decrease. During some time, 
the two surfaces will have a common point (Fig. 1). 


Sik (x, Ry, t) = 0 
Xi 


SA 


T (x, Ry, t) =0 


AG 


Fig. 1. Geometric interpretation of hypersurfaces in phase space under guidance in boost phase 


At the point S; and 7°, they have a common tangency, that is, their gradients coincide in the phase space. In this 
way: 

- existence of the optimal trajectory x,,,(¢) and optimal control u,,, (¢) is noted; 

- termination of the boost phase of the trajectory (¢ = ¢,) and the transition to the ballistic flight phase are deter- 
mined. 

At the tangency point, the isochronous surface and the hypersurface of the boundary conditions have a com- 
mon tangency and a normal. Mathematically, the condition for the existence of a common normal is determined by the 


expression [2]: 
0 
or") _ 4{@.\. a 
Ox, Ox, 


where A is constant determined from the analysis of the convexity of both hypersurfaces. 
From the relations (11) and (18), we transform (17). We obtain the expression for the direction cosines of the 
thrust vector as a function of the current and final motion parameters: 
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4- (—2x,x,%, eee 2x; X, oe ) (19) 


2 2 9 
2 2 2 1/2 
(2x, x3, Nay Nag 20 —x{X,, ) + (2x,x LN Ky Nyt LN XN Py Ny Nyy +X,%;%,, } | 


2 
(2x; Ky 2X NX, PLN HX NHK, FX NX, ) 


B- (20) 


2 2 1/2 ? 
2 2 2 
| (2x, x,3, +X,X,X,, +2X,X, —X4%,, ) + (2x,%, LN kh 2K hin NN, +X,%5%,, } | 


where 4=cosa,, B=cosa,. 
From the relations (19), (20), it is easy to obtain a parameter natural for this UAV type to determine the orien- 


tation vector of the control acceleration (in a certain sense, of the thrust vector), the pitching angle: 
COS, _ Coes = 20K N, 2H X,, PX AN, FH MX, ) (21) 


O(S) = arctg 
(-2x,x,%, to tO a, ) 





O 
COS QL, 

In our case, the UAV motion takes place in the given plane, and the yaw angle is  (S') = 0. Using the proposed 
methodology for the synthesis of terminal optimal control, an algorithm for calculating the pitching angle was devel- 
oped [2], and the computational simulation of the UAV flight was carried out [7]. 

Research Results. The numerical studies have been conducted using the software that implements the pro- 
posed method. We are talking about the algorithmic software for terminal guidance of ballistic aircraft based on the 
solution to boundary problems of ballistics. The corresponding computer program was registered in 2013. 

When modeling, a hypothetical accelerating tool was used with the characteristics and initial conditions for the 


UAV launch ascent, given in [8]. 
The initial conditions of the simulation and the calculation results of the early ascent are given in Table 1. 


Table 1 
Initial conditions for UAV launch ascent 


xin 
5.6323 —0.0000 322.6757 90.00 0.00 0.000 
5.6323 11.1626 322.6756 90.00 1.98 —0).048 
21.4565 69.4705 322.6727 69.34 12.13 —7.881 
155.6274 259.0827 322.6335 60.00 35.9] 0.011 


368.3725 464.9493 322.544] 50.02 43.59 —2.076 
662.2567 679.2317 322.4053 43.44 42.79 —2.582 
1044.7345 915.1428 322.2176 38.74 40.11 —2.677 
1429.1498 1131.2610 322.0334 36.01 37.94 —2.501 





Fig. 2 shows the pitching-angle profile in the boost flight phase. 
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Dependence of pitching and attack angles on apparent speed rate 





0 1 000 2 000 3000 4000 3000 6000 F000 


— pitching angle — attack angle 


Dependence of UAV entrance angle on apparent speed rate 





0 1 000 2 000 3000 4000 3 000 6 000 F000 


— UAV entrance angle 
Fig. 2. Variation of UAV pitching and entrance angles 


Fuller information on the simulation results given in [8] is as follows. 

1. The resulting control is suboptimal due to the use of boundary conditions in an analytical form. To improve 
the accuracy of the boundary conditions in the navigation algorithm, it 1s necessary to introduce an up-date algorithm, 
which is based on attracting more accurate models of the Earth’s gravitational field. Hence, the launch ascent accuracy 
will be improved. 

2. The effect of random disturbances is compensated by the adaptive properties of terminal guidance, as well 
as by predicting the motion parameters and the formation of control at a decreasing time interval each time until the end 
of the launch ascent process. 

3. The obtained estimates under modeling allowed us to rationally select the discreteness of the navigation 
and guidance algorithms and thereby limit the range of requirements for the FMC during its implementation. 

4. The implementability of guidance algorithms based on the developed methodology on modern onboard 
computers creates no difficulties. The required response rate is (1-1.5)x10° k.o./s. 

Discussion and Conclusions. Thus, the principle of dynamic programming seems to be most preferable. It 
should be used when solving the problem on the centroidal motion control synthesis (guidance problem) of the UAV 
with long-range capabilities in the boost phase. The well-known remark on the applicability of dynamic programming, 
the so-called “curse of dimensionality”, is inappropriate in the task of developing control as a function of the final state 
[5]. In addition, the use of boundary conditions in an analytical form simplifies essentially the formation of suboptimal 
control and enables to change strategically the flight missions. This scales up the applicability of this algorithm for 
UAV of various purposes. 
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